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ABSTRACT 


The  primary  goal  of  this  projeet  is  to  improve  the  eapability  to  identify  and  deteet  surfaee  waves  for  the 
purpose  of  earthquake/explosion  diserimination.  We  are  developing  improved,  higher  resolution  earth  and 
dispersion  models.  The  models  eonsist  of  approximately  600  distinet  erust  and  upper  mantle  struetures, 
with  surfaee  layering  and/or  oeean  depths  that  vary  on  a  one-degree  grid.  There  are  a  total  of  64,800  earth 
models  and  dispersion  eurves,  but  the  tomographie  inversion  is  performed  only  for  the  600  distinet  erust 
and  upper  mantle  models,  with  the  shallow  strueture  eonstrained  by  other  information.  The  data  set  used  in 
the  inversion  now  eonsists  of  approximately  548,000  phase  and  group  veloeity  dispersion  measurements 
obtained  from  a  variety  of  sourees.  Surfaee  sediments  are  defined  using  the  global  sediment  maps  of  Laske 
and  Masters  (1997),  and  oeean  bathymetry  is  defined  using  the  Etopo5  topographie  data  set. 

Automatie  identifieation  of  surfaee  waves  at  the  International  Data  Centre  is  eurrently  performed  by 
narrow-band  filtering  the  data  at  several  frequeneies,  and  then  eomparing  the  arrival  times  with  a 
regionalized  dispersion  model.  We  have  implemented  and  tested  a  new  proeedure  in  whieh  we  first  phase- 
mateh  filter  the  data  and  then  apply  narrow-band  filters  to  the  eompressed  waveform  and  use  a  deteetion 
test  similar  to  the  eurrent  test.  This  allows  us  to  take  advantage  of  the  improved  signal-to-noise  ratio  of  the 
phase-mateh  filtered  waveforms,  while  retaining  the  robustness  of  narrow-band  filtering  for  frequeney- 
dependent  signal  identifieation.  After  phase-matehed  filtering,  the  predieted  arrival  time  is  zero  at  all 
frequeneies,  so  we  test  to  see  if  the  arrivals  are  within  a  time  window  similar  to  that  used  in  the  existing 
test.  To  test  the  proeedure,  we  proeessed  the  same  data  set  using  five-degree  and  one-degree  models  with 
and  without  phase-matehed  filtering.  Deteetions  using  the  one-degree  model  with  phase-matehed  filtering 
inereased  by  more  than  30%  eompared  to  the  five-degree  model  eurrently  in  use  at  the  IDC  without  phase- 
matehed  filtering. 

We  use  long-period  waveforms  from  historie  nuelear  explosions  to  assess  the  potential  of  automated 
Rayleigh  wave  travel-time  pieks  to  improve  seismie  event  loeations.  The  improved  aeeuraey  of  loeations 
reported  by  Yaeoub  (2000),  based  on  20-seeond  Rayleigh  waves,  provided  the  impetus  for  this  work. 
Although  surfaee  wave  arrival  times  eannot  be  measured  as  aeeurately  as  body  wave  arrivals,  that  surfaee 
waves  are  mueh  slower  means  that  aeeurate  loeations  ean  be  aehieved.  We  find  that  good  loeations  ean  be 
determined  with  surfaee  waves,  whieh  eould  potentially  improve  loeations  made  with  body  waves  alone 
provided  that  the  paths  are  relatively  short  and  an  aeeurate  dispersion  model  is  available  for  the  region. 

This  may  be  espeeially  important  for  small  events  with  few  total  measurements.  To  assess  the  value  of 
single  surfaee  wave  measurements  for  that  purpose,  we  also  estimate  the  aeeuraey  of  single -station  distanee 
estimates. 

KEY  WORDS:  surfaee  wave,  dispersion  eurve,  phase-matehed  filter,  regionalization,  moment,  loeation 


OBJECTIVE 

The  primary  objeetive  of  this  projeet  is  to  develop  improved  methods  for  deteetion  and  measurement  of 
surfaee  waves.  This  is  important  beeause  the  threshold  level  for  surfaee  waves  is,  in  general,  the  threshold 
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for  the  Msirrib  discriminant.  By  reducing  the  threshold  for  surface  wave  detection,  we  can  reduce  the 
number  of  unidentified  events  and  increase  the  likelihood  of  identifying  smaller  nuclear  explosions. 

RESEARCH  ACCOMPLISHED 


As  stated  above,  the  primary  focus  of  this  project  is  to  improve  the  capability  to  identify  and  detect  surface 
waves  for  the  purpose  of  earthquake/explosion  discrimination.  This  is  being  accomplished  through 
development  of  high-resolution  earth  models  constrained  by  inversion  of  surface  wave  observations  and 
auxiliary  information,  which  in  turn  are  used  to  derive  surface  wave  phase  and  group  velocity  dispersion 
curves  on  a  global  grid.  These  dispersion  curves  are  then  used  to  construct  phase-matched  filters  to  improve 
surface  wave  detection.  The  availability  of  high-resolution  group  velocity  curves  also  makes  it  possible  to 
use  surface  waves  to  improve  event  location,  as  suggested  by  Yacoub  (2000),  so  we  have  performed  a 
study  to  assess  under  what  conditions  this  technique  may  be  useful.  This  paper  is  organized  into  three  parts: 
improvement  of  regionalized  global  earth  models  and  dispersion  curves;  automatic  surface  wave 
identification  using  phase-matched  filtering;  and  location  of  events  using  surface  wave  observations. 

Development  of  Regionalized  Earth  and  Dispersion  Models 

The  most  important  information  required  for  improving  surface  wave  detection  and  measurement  is 
accurate  global  dispersion  maps.  Consequently,  most  of  our  work  has  concentrated  on  the  development  of 
three-dimensional  velocity  models  of  the  earth’s  crust  and  upper  mantle  by  tomographic  inversion  of 
surface  wave  dispersion  measurements.  This  work  uses  an  extension  of  the  technique  described  by  Stevens 
and  McLaughlin  (2001),  in  which  a  global  earth  model  with  149  distinct  model  types  on  a  five-degree  grid 
was  developed.  That  model  is  now  being  used  for  routine  surface  wave  identification  at  the  International 
Data  Centre  (IDC).  The  technique  used  at  the  IDC  is  to  compare  predicted  group  velocity  dispersion  curves 
derived  from  these  regionalized  models  with  measured  dispersion  curves  from  observed  surface  waves.  As 
discussed  later,  phase  velocities  derived  from  these  earth  models  can  also  be  used  to  develop  phase- 
matched  filters  to  improve  signal-to-noise  ratio  and  optimize  the  detection  test. 

The  data  used  to  derive  our  earth  models  consist  of  more  than  548,000  measurements  of  surface  wave 
group  and  phase  velocity  dispersion  from  earthquakes  and  explosions.  These  measurements  were  obtained 
from  various  sources  and  are  for  frequencies  ranging  between  0.005  and  0.1667  Hz  with  the  great  majority 
of  observations  between  0.01  and  0.1  Hz.  The  data  set  has  been  derived  from  a  variety  of  regional  and 
global  studies  including  the  following:  global  surface  wave  group  velocities  from  earthquakes  derived 
using  PIDC  GSETT3  data  (Stevens  and  McLaughlin,  1996),  augmented  with  more  recent  measurements 
derived  from  PIDC  data;  surface  wave  phase  and  group  velocity  dispersion  curves  from  underground 
nuclear  test  sites  (Stevens,  1986;  Stevens  and  McLaughlin,  1988),  calculated  from  earth  models  for  270 
paths  (test  site  -  station  combinations)  at  10  frequencies  between  0.015  and  0.06  Hz;  phase  and  group 
velocity  measurements  for  western  Asia  and  Saudi  Arabia  from  Mitchell  et  al.(1996)  for  12  paths  at  17 
frequencies  between  0.012  and  0.14  Hz;  the  global  phase  velocity  model  of  Ekstrom  et  al.  (1996)  for  9 
periods  between  35  and  150  seconds  calculated  for  each  5-degree  grid  block  from  a  spherical  harmonic 
expansion  of  order  1=40;  group  velocity  measurements  for  Eurasia  from  Ritzwoller  et  al.(1996)  and 
Levshin  et  al.(1996)  for  20  frequencies  between  0.004  and  0.1  Hz  with  500  to  5000  paths  per  frequency; 
Antarctic  and  South  American  group  velocity  measurements  from  the  University  of  Colorado  (Vdovin  et 
al.,  1999;  Ritzwoller  et  al.,  1999);  high-frequency  Eurasian  dispersion  measurements  from  University  of 
Colorado  (Levshin  and  Ritzwoller,  pers.  comm.,  2001),  and  a  large  set  of  dispersion  measurements  from 
Saudi  Arabia  provided  by  Herrmann  and  Mokhtar  at  St.  Louis  University.  We  have  received  a  large  set  of 
data  from  LLNL  for  North  Africa  and  Eurasia  (Pasyanos,  pers.  comm.,  2001),  but  are  still  reviewing  the 
measurements  and  have  not  included  it  in  the  inversion  results  described  here. 

Development  of  improved  earth  models  and  dispersion  curves  has  proceeded  using  the  following  approach: 

1.  The  starting  point  was  the  5 -degree  IDC  149  model  set  (Stevens  and  McLaughlin,  2001),  which  was 
based  on  approximately  90,000  dispersion  measurements. 

2.  New  dispersion  measurements  were  added  to  the  data  set  and  the  same  model  set  was  reinverted  with 
the  new  data. 
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3.  New  model  types  were  added  in  areas  with  inereased  data  or  where  the  data  misfit  indieated  that  new 
model  types  were  required  and  a  new  inversion  was  performed  with  the  larger  data  and  model  set. 

4.  A  proeedure  was  developed  for  ineluding  shallow  strueture,  partieularly  sediment  thieknesses  and 
oeean  depths,  on  a  one  degree  grid  while  inverting  for  a  limited  set  of  distinet  models,  most  of  whieh 
remained  on  a  five-degree  grid,  in  the  erust  and  upper  mantle. 

5.  The  loeations  and  boundaries  of  the  underlying  model  types  were  redefined,  starting  with  the  Crust  2.0 
model,  so  that  they  follow  plate  boundaries  and  other  geologie  eonstraints.  The  inversion  teehnique 
was  also  modified  to  allow  small  variations  in  Moho  depth  and  layer  thiekness  to  be  defined  on  a  one- 
degree  grid  for  eaeh  model  type. 

So,  with  our  eurrent  approaeh,  the  inversion  is  performed  for  shear  veloeity  struetures  in  approximately  500 
distinet  erust  and  upper  mantle  model  types,  but  the  shallow  strueture,  bathymetry,  as  well  as  small  ehanges 
in  layer  thiekness  and  Moho  depth,  ean  vary  on  a  one-degree  grid  (64,800  distinet  models).  P  wave 
veloeities  are  eonstrained  via  a  eonstant  Poisson’s  ratio  of  0.27,  and  density  via  Bireh’s  Law.  The  layers 
vary  in  thiekness  (averaging  about  10  km)  and  reaeh  down  to  a  depth  of  about  250  km.  Below  250  km  the 
Earth  model  is  fixed  to  mateh  PREM.  The  top  few  kilometers  of  the  model  (eonsisting  of  water,  iee  and/or 
sediments)  is  fixed  and  matehes  data  from  1 -degree  bathymetry  maps  made  by  averaging  Etopo5  5  minute 
measurements  of  topography,  and  Laske  and  Masters  (1997)  1 -degree  maps  of  sediments.  One  advantage 
of  inverting  for  a  fixed  set  of  erust  and  upper  mantle  struetures  is  that  it  reduees  the  problem  to  a  more 
manageable  size.  An  earth  model  eonsisting  of  16200  2-by-2  degree  eells  and  with  about  10  layers  per  eell 
would  have  on  the  order  of  10^  free  parameters.  Our  eurrent  model  eonsists  of  537  model  types  thatreduee 
the  total  number  of  free  parameters  to  6494.  Other  advantages  of  this  type  of  inversion  are  that  it  uses  all 
frequeneies  (and  both  phase  and  group  veloeity)  at  onee,  with  different  frequeneies  resolving  different 
length  seales;  the  teehnique  keeps  similar  struetures  eonsistent  with  eaeh  other,  avoiding  fluetuations  in 
nearby  points  eommon  in  more  traditional  group  veloeity  inversions;  and  it  allows  Moho  depths,  suberustal 
layer  thieknesses,  and  lateral  boundaries  between  geologie  regions  to  be  added  to  the  model  as  fixed  a 
priori  information.  The  prineipal  disadvantages  are  that  distinet  model  types  must  be  ehosen  earefully  and 
added  sparingly  if  needed,  and  that  there  is  no  smoothness  eonstraint  between  adjaeent  struetures. 

The  eurrent  model  under  development  eonsists  of  537  model  types.  There  are  two  main  varieties  of  types, 
those  originating  from  the  Crust  2.0  2x2  degree  erustal  types  (Bassin  et  al.,  2000  and  Laske  et  al.  2001)  and 
those  based  on  oeean  ages  (Stevens  and  Adams,  2000).  The  Crust  2.0  models  are  an  improvement  of  the 
Mooney  et  al.  (1998)  Crust  5.1  model  revised  and  refined  to  a  2-degree  seale.  They  may  inelude  a  layer  of 
iee,  water,  1-3  layers  of  sediments,  and  3  erustal  layers.  Beeause  many  of  the  Crust  2.0  models  differ  only 
in  shallow  strueture  rather  than  erustal  veloeity  or  thiekness,  we  use  a  somewhat  different  set  of  models. 

We  eombine  similar  model  types  that  vary  in  shallow  strueture,  sinee  we  treat  those  variations  as 
independent  eonstraints,  whieh  has  the  effeet  of  redueing  the  number  of  model  types  relative  to  Crust  2.0. 
However,  we  also  separate  models  in  widely  separated  regions  so  that  results  from  a  model  in  North 
Ameriea,  for  example,  do  not  affeet  a  model  of  the  same  type  in  Eurasia.  This  inereases  the  number  of 
model  types.  We  also  reparameterized  the  strueture  of  the  oeeans  based  on  oeean  age.  Previously  oeean 
struetures  were  based  on  both  age  and  sedimentary  thiekness.  We  formed  17  oeean  types  based  on  the 
Muller  et  al.  (1997)  isoehron  map.  There  are  five  types  for  eaeh  of  the  three  major  oeeans:  Paeifie,  Indian 
and  Atlantie.  The  five  types  mateh  the  age  ranges  0  to  10.9  Ma,  10.9  to  20.1  Ma,  20.1  to  40.1  Ma,  40.1  to 
83.5  Ma  and  83.5  to  180.0  Ma.  The  remaining  two  types  are  for  the  Aretie  Sea  for  ages  0  to  40.1  Ma,  and 
40.1  to  180.0  Ma. 

Sedimentary  struetures,  bathymetry  and  lateral  boundaries  between  different  oeeanie  types  are  ineluded  on 
a  1 -degree  level  of  detail.  This  parameterization  serves  as  an  initial  mo  del  for  a  series  of  3-D  tomographie 
inversions.  Regularization  is  eontrolled  by  two  parameters,  one  for  a  vertieal  smoothing  eondition  applied 
to  layers  in  eaeh  model  type,  and  a  damping  parameter.  The  tomographie  inversion  has  the  form  shown  in 
equation  1 .  Sg  is  the  slowness  of  the  g^^  one  degree  eell;  mj  is  the  shear  wave  veloeity  of  the  layer  where) 
spans  all  the  adjustable  layers  of  all  model  types.  Xig  is  the  length  of  i^^  raypath  traveled  in  the  g^^  eell,  ?  Ti 
is  the  travel  time  residual  for  the  i^^  raypath.  From  this  equation  we  ean  see  that  64800  one-degree  eells  are 
needed  to  ealeulate  the  partial  derivatives.  One-degree  eells  are  neeessary  here  to  earry  the  information 
about  sediments  and  bathymetry.  However,  the  free  parameters,  indexed  with),  span  a  eoarser  seale,  that  of 
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the  model  types.  These  equations,  plus  equations  for  dampening  and  smoothing,  are  solved  using  the  LSQR 
method.  20  send  phase  and  group  veloeity  maps  are  shown  in  figures  1  and  2. 


Figure  1.  20-seeond  group  veloeity  map. 


3.0  3^  3.4  3.0  3.3  4.D 


Figure  2.  20-seeond  phase  veloeity  map 


The  depth  profiles  of  the  three  main  oeeanie  types  are  shown  in  figure  3 . 
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Figure  3.  Shear  wave  veloeity  depth  profiles  for  the  main  fifteen  oeean  types.  These  results  are  eonsistent 
with  thermal  models  of  oeeanie  lithosphere. 


Automatic  Surface  Wave  Detection  Using  Phase-Matched  Filtering 

Automatic  identification  of  surface  waves  at  the  International  Data  Centre  is  currently  performed  using  the 
processing  program  Maxsurf  by  narrow-band  filtering  the  data  at  several  frequencies  and  then  comparing 
the  arrival  times  with  a  regionalized  dispersion  model.  An  automatic  surface  wave  processing  program, 
Maxpmf,  has  been  developed  which  is  similar  to  Maxsurf  except  that  it  applies  phase-matched  filtering  to 
seismograms  and  calculates  path-corrected  spectral  magnitudes  in  addition  to  Mg.  Maxpmf  integrates  a 
regionalized  phase  velocity  model  to  generate  a  phase-matched  filter  that  is  used  to  compress  the  surface 
wave  waveforms.  Figure  4  below  shows  an  example  of  waveform  compression  after  application  of  phase- 
matched  filters  to  a  data  set.  Figure  5  shows  the  recommended  way  to  perform  surface  wave  identification 
using  phase-matched  filtering,  which  is  to  apply  narrow-band  filters  to  the  phase-match  filtered  waveforms 
and  then  look  for  arrival  times  near  zero. 
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Figure  4.  Data  in  the  5  km/sec  to  2  km/sec  time  window  from  an  mb  3.9  South  Pacific  earthquake  after 

transformation  to  a  KS36000  instrument  (left)  and  after  compression  by  phase-matched  filtering 
(right). 
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Figure  5.  The  narrow-band  filter  deteetion  test  (top)  eonpares  measured  group  veloeity  with  a  regionalized 
group  veloeity  model.  To  improve  deteetion,  we  phase-mateh  filter  the  data  first,  then  apply 
narrow-band  filters  (bottom).  The  deteetion  test  is  applied  for  the  same  time  interval,  but  the  time 
is  now  eentered  around  zero.  This  example  is  for  an  mb  4.2  South  Ameriean  earthquake  observed 
at  BDFB,  BOSA,  and  ULM. 

As  a  test  of  both  the  improvement  in  models  and  the  ability  of  phase-mateh  filtering  to  improve  deteetion, 
we  ran  Maxpmf  on  one  full  day  of  data  and  applied  both  the  eurrent  narrow-band  filtered  test  and  the 
phase-matehed  filter  test  deseribed  above  to  this  data  set.  The  results  are  shown  in  Figure  6.  The  four  bars 
show,  respeetively,  the  number  of  deteetions  determined  by  the  eurrent  proeedure,  the  nurrber  of  deteetions 
with  the  best  5 -degree  model,  the  number  of  deteetions  from  the  eurrent  one-degree  model,  and  the  nurrber 
of  deteetions  using  the  one-degree  model  with  phase-matehed  filtering.  Eaeh  improved  model,  and  the 
implementation  of  phase-matehed  filtering,  shows  an  improvement.  The  number  of  deteetions  inereased  by 
32%  over  eurrent  IDC  deteetions  with  phase-matehed  filtering  using  the  most  reeent  one-degree  model. 

The  one-degree  model  used  in  this  test  was  a  modifieation  of  the  earlier  5 -degree  model  and  did  not  inelude 
the  modifieations  diseussed  earlier  to  ineorporate  the  Crust  2.0  improvements. 
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Figure  6.  Improvement  in  the  number  of  deteetions  for  a  day  of  data  with  improvements  in  the  model,  and 
with  phase-matehed  filtering. 
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As  an  additional  test  of  the  models,  we  eompare  observed  and  predieted  dispersion  eurves  for  data  from 
nuelear  test  sites.  Tables  1  and  2  show  group  veloeity  dispersion  residuals  and  standard  deviations  in  the 
residuals  for  periods  of  40  seeonds  and  20  seeonds,  respeetively,  for  paths  from  several  nuelear  test  sites. 
With  the  expanded  data  set  and  149  models,  the  results  are  ehanged  only  slightly  relative  to  the  inversion 
with  the  smaller  data  set.  In  faet  the  results  are  degraded  slightly  for  the  paths  from  the  test  sites.  This  is 
beeause  the  new  data  eaused  larger  ehanges  in  other  areas  and  with  the  eonstraint  of  the  small  number  of 
models,  the  data  fit  beeame  slightly  worse  in  order  to  allow  larger  improvements  in  the  areas  with 
eoneentrations  of  new  data.  With  388  models,  however,  these  eonstraints  are  relaxed,  and  the  data  fit 
improves  for  all  test  sites,  and  improves  quite  dramatieally  for  paths  in  Asia.  Curiously,  the  data  fits  are  not 
improved  with  the  one-degree  model.  This  result  together  with  the  result  shown  in  Figure  8  suggests  that 
the  global  model  is  improved,  but  at  the  eost  of  some  degradation  in  paths  from  the  major  test  sites.  We  are 
in  the  proeess  of  investigating  the  reasons  for  this  now. 


Table  1.40  Seeond  Group  Veloeity  %  Average  Residuals  (Standard  Deviations) 


Source 

IDC5° 

149  Models 

5°  149  Models 
Expanded  Data  Set 

5° 

388  Models 

1° 

489  Models 

NTS  (59) 

East  Kazakh  (40) 
Mururoa  (13) 

Novaya  Zemlya  (99) 
Amchitka  (55) 

0.15(1.40) 

0.80  (2.21) 
-0.78  (1.43) 
-0.06  (2.06) 

0.25  (1.39) 

0.22(1.44) 

0.85  (2.45) 

-0.73  (1.79) 

0.01  (2.05) 

0.34  (1.39) 

0.14(1.20) 

0.13  (1.11) 

-0.59  (1.56) 
-0.18(1.81) 

0.35  (1.08) 

-0.17(1.28) 

0.31  (1.18) 
-1.07(1.48) 
-0.10(1.98) 
0.08(1.18) 

Table  2.20  Seeond  Group  Veloeity  %  Average  Residuals  (Standard  Deviations) 


Source 

IDC5° 

149  Models 

5°  149  Models  Expanded 
Data  Set 

5° 

388  Models 

1° 

489  Models 

NTS  (58) 

-0.47  (2.25) 

-0.38  (2.23) 

-0.30  (1.88) 

-0.08  (2.61) 

East  Kazakh  (40) 

-0.75  (1.96) 

-0.65  (2.07) 

-0.64(1.44) 

-0.05  (1.27) 

Mururoa  (11) 

0.41  (1.43) 

0.41  (1.40) 

0.53  (1.36) 

-0.05  (2.22) 

Novaya  Zemlya  (99) 

-0.26  (3.61) 

-0.23  (3.66) 

-0.43  (3.01) 

0.31  (3.28) 

Amchitka  (54) 

0.83  (3.76) 

0.70  (3.74) 

0.22  (3.23) 

0.84  (3.92) 

Location  Improvement  Using  Surface  Waves 

Surface  wave  arrival  times  can  be  used  in  the  same  manner  as  (or  together  with)  P-wave  arrival  times  to 
determine  source  location.  With  a  good  regionalized  group  velocity  model  and  well-determined  group 
arrival  times,  the  location  could  be  determined  quite  accurately,  and  in  cases  where  P-wave  coverage  is 
poor,  or  azimuthal  coverage  is  limited,  using  surface  wave  arrivals  has  the  potential  to  significantly 
improve  location  accuracy.  Yacoub  (2000)  found  that  he  was  able  to  determine  the  location  of  7  Nevada 
Test  Site  (NTS)  explosions  more  accurately  with  surface  waves  narrow-band  filtered  at  20  seconds  than 
with  P  waves.  This  was  done  using  a  constant  group  velocity  of  3.0  km/sec  for  all  paths.  This  result  is 
attributed  to  the  fact  that  surface  wave  group  velocities  are  much  slower  than  P-wave  velocities,  so  that 
errors  in  arrival  time  correspond  to  significantly  smaller  errors  in  distance  than  the  corresponding  errors  in 
P  wave  arrival  times,  and  to  the  consistency  in  surface  wave  arrival  times. 

We  have  performed  some  experiments  using  our  global  group  velocity  models  and  dispersion 
measurements  to  assess  the  ability  to  locate  events  using  surface  wave  arrivals.  Table  3  lists  location  errors 
from  Yacoub  (2000)  derived  using  both  surface  waves  and  P  waves,  and  the  location  errors  we  derived 
from  independent  measurements  of  data  at  several  frequencies  from  the  same  events  using  the  regionalized 
dispersion  curves  discussed  earlier.  The  results  confirm  that  location  estimates  are  somewhat  better  than 
can  be  obtained  from  P-wave  arrivals  at  the  same  stations,  although  the  location  errors  are  somewhat  larger 
than  obtained  by  Yacoub.  This  may  be  because  Yacoub  used  a  larger  and  higher  quality  data  set.  Also,  it 
should  be  recognized  that  Yacoub ’s  P-wave  measurements  used  only  data  from  stations  that  also  had 
surface  wave  data,  and  considerably  better  results  could  be  obtained  from  a  larger  P-wave  data  set. 

Yacoub  (2000)  found  remarkably  constant  group  velocity  measurements,  all  very  close  to  3.0  km/sec.  We 
cannot  confirm  this  and  assume  that  it  must  be  due  either  to  a  particular  choice  of  paths,  or  some  anomaly 
in  the  way  the  measurements  were  made.  Figure  7  shows  the  distribution  of  measured  group  velocities  from 
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all  of  the  explosion  data  from  all  test  sites.  While  the  mean  of  the  distribution  is  very  elose  to  3.0  km/see, 
there  is  eonsiderable  variation  from  this  value  in  the  data  set. 


Table  3.  Loeation  errors  for  3  methods,  the  number  of  stations  used  for  eaeh  (Yaeoub  P  and  surfaee  wave 
loeations  used  identieal  sets  of  stations). 


Event 

#STA 

Loeation  Error  (km) 

Yaeoub 

IDC 

Yae.  LR 

Yae.P 

IDC 

Hearts 

19 

17 

11.2 

18.3 

22.5 

Baekbeaeh 

19 

2 

9.7 

33.1 

NA 

Seantling 

18 

11 

25.1 

49.9 

18.9 

Lowball 

17 

12 

16.3 

28.3 

18.1 

Mizzen 

14 

9 

9.8 

32.0 

34.2 

Strake 

14 

0 

13.0 

26.0 

NA 

Sheepshead 

12 

16 

9.2 

16.1 

13.0 

Averages 

13  ±6 

29+  11 

21  ±7 

m 


Sronj  Vdocity 

Figure  7.  Histogram  of  measured  group  veloeities  of  20-seeond  Rayleigh  waves  from  explosions  at  all 
major  test  sites. 
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Figure  8.  Mean  and  two  standard  deviation  error  bounds  of  distanee  estimates  from  single  stations,  for 
different  frequeneies.  There  are  approximately  100  data  per  distanee  bin. 
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Figure  8  shows  the  estimated  error  as  a  funetion  of  distanee  for  single  station  measurements.  In  general, 
the  loeation  error  inereases  with  distanee  to  the  observing  station.  Lower  frequeneies,  however,  degrade 
mueh  more  slowly.  However,  at  shorter  distanees  higher  frequeney  data  give  more  aeeurate  loeations. 

Figure  9  shows  the  loeation  error  derived  from  surfaee  wave  measurements  for  the  Semipalatinsk  test  site 
as  a  funetion  of  the  number  of  reeording  stations.  The  results  are  similar  to  the  NTS  results  with  a  typieal 
loeation  error  of  about  30  km.  It  is  important  to  note,  however,  that  loeation  aeeuraey  depends  as  mueh  or 
more  on  whieh  stations  reeorded  the  event  as  on  how  many.  Referring  to  Figure  8,  it  is  elear  that  surfaee 
waves  reeorded  at  distanees  of  3000  km  or  less  provide  mueh  better  eonstraints  on  loeation  than  more 
distant  stations.  Surfaee  waves  are  most  likely  to  be  useful  in  eonstraining  loeation  therefore  when  they  are 
measured  at  shorter  distanees  and  eombined  with  a  small  number  of  P-wave  measurements.  The  P-wave 
measurements  on  their  own  may  not  be  very  aeeurate,  but  even  a  single  surfaee  wave  aeeurately  measured 
at  a  distanee  of  <2000  km  eould  signifieantly  improve  the  loeation. 


Figure  9.  Loeation  errors  for  Semipalatinsk  explosions.  The  median  error  for  events  reeorded  at  6  or  more 
stations  is  30  km. 

CONCLUSIONS  AND  RECOMMENDATIONS 

This  study  has  foeused  on  development  of  improved  methods  for  deteetion  and  measurements  of  surfaee 
waves.  As  part  of  this  effort,  we  have  developed  a  gradually  improving  set  of  global  earth  models  suitable 
for  ealeulation  of  phase  and  group  veloeity  dispersion  eurves.  Over  the  eourse  of  this  study,  the  earth 
models  have  evolved  from  a  five-degree  model  based  on  90,000  dispersion  measurements  to  a  one-degree 
model  based  on  over  500,000  measurements.  Currently,  we  are  ineorporating  improvements  in  boundaries 
and  eonstraints  on  shallow  strueture  and  Moho  depth  using  the  Crust  2.0  models  to  improve  the  starting 
models  for  the  tomographie  inversions.  We  have  used  phase  veloeities  derived  from  these  models  to  form 
phase-matehed  filters,  and  developed  an  improvement  on  the  surfaee  wave  deteetion  algorithm  that  uses 
phase-matehed  filtering  followed  by  narrow-band  filtering  for  deteetion.  We  performed  a  test  on  one  day  of 
data  at  the  PIDC  with  several  dispersion  models  and  found  a  signifieant  improvement  using  this  teehnique. 
We  have  also  tested  the  use  of  surfaee  waves  for  loeation.  Although  we  eannot  eonfirm  improvement  as 
signifieant  as  suggested  by  Yaeoub  (2000),  it  does  appear  that  surfaee  waves  ean  add  an  additional 
eonstraint  that  eould  signifieantly  improve  loeation,  partieularly  in  eases  where  there  are  only  a  few  P 
arrivals  and  the  surfaee  waves  are  measured  on  short  paths. 
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